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Abstract. Semi-infinite d-dimensional systems with an m-axial bulk Lifshitz 
point are considered whose (d — l)-dimensional surface hyper-plane is oriented 
perpendicular to one of the m modulation axes. An n-component (j> 4 field theory 
describing the bulk and boundary critical behaviour when (i) the Hamiltonian 
can be taken to have 0(n) symmetry and (ii) spatial anisotropics breaking 
its Euclidean symmetry in the m-dimensional coordinate subspace of potential 
modulation directions may be ignored is investigated. The long-distance 
behaviour at the ordinary surface transition is mapped onto a field theory 
with the boundary conditions that both the order parameter <p and its normal 
derivative 9 n vanish at the surface plane. The boundary-operator expansion 
is utilized to study the short-distance behaviour of <p near the surface. Its 
leading contribution is found to be controlled by the boundary operator c9 2 <fi. The 
field theory is renormalized for dimensions d below the upper critical dimension 
d*(m) = 4 + m/2, with a corresponding surface source term oc d%<fi added. The 
anomalous dimension of this boundary operator is computed to first order in 
e = d* — d. The result is used in conjunction with scaling laws to estimate the 
value of the single independent surface critical exponent P^ d,± '^ for d = 3. Our 
estimate for the case m = n = 1 of a uniaxial Lifshitz point in Ising systems is in 
reasonable agreement with published Monte Carlo results. 

E-mail: shpotaph.icmp.lviv.ua 

1. Introduction 

The significance of the n-component 4 models with 0(n) symmetric Hamiltonian 
derives from the fact that they represent the most common — and probably also 
most important — universality classes of critical behaviour at bulk critical points of 
d-dimensional systems with short-range interactions. Prominent examples of such 
universality classes for given d and n = 1, 2, and 3 are those of the Ising, XY, and 
isotropic Heisenberg models, respectively. 

When such systems are bounded by free (d — l)-dimensional hyper-surfaces or 
walls, a wealth of distinct boundary critical phenomena can occur P]|21E1E]- A well- 
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studied case are systems with a free surface that can be described by a semi-infinite 
n-component </> 4 model. The Hamiltonian of the latter is of the form 

H= [ C h (x)dV+ [ A(r)dA, (1) 

where and J B mean volume and surface integrals over the half-space 2J = = 
{x = (r,z)\r G R d_1 ,0 < z < 00} and the z = boundary hyperplane 23, 
respectively. Provided neither bulk nor surface terms breaking the 0(n) symmetry 
must be included, the bulk and surface densities are given by 

A = ^(V0) 2 + ^ 2 + ||0| 4 (2) 

and 

A = ^ 2 . (3) 

Furthermore, one can distinguish three different types of surface transitions that take 
place at the bulk critical temperature T c They are called ordinary, special, and 

extraordinary surface transitions and occur, depending on whether c is larger than, 
equal to, or smaller than a certain critical value c sp . 

The extraordinary transition is a transition from a surface-ordered, bulk- 
disordered high-temperature phase to a low-temperature phase with long-range order 
both at the surface and in the bulk. It — and hence also the special transition — can 
occur only for those choices of d and n for which fluctuations do not destroy long-range 
order at the surface at all nonzero temperatures T . This requires d — 1 > 1 in the 
scalar case n = 1 , and d — 1 > 2 in the continuous-symmetry case n > 1. 

Suppose n and d are such that all three of these transitions can occur. Then the 
critical behaviour that surface quantities exhibit at any of these is representative of a 
separate surface universality class, although their bulk critical behaviour is the same. 
Thus the bulk universality class associated with this choice of n and d splits into three 
distinct surface universality classes. 

In the present paper we shall be concerned with the surface critical behaviour of n- 
vector systems at m-axial bulk Lifshitz points. A Lifshitz point (LP) is a multicritical 
point at which a disordered, a homogeneous ordered, and a modulated ordered phase 
meet EH H| ■ A family of natural extensions of the bulk models defined on M d by the 
bulk density J5J) that have such LPs was introduced decades ago [El El, but investigated 
in greater detail via field-theoretic renormalization group (RG) only in the past few 
years (see the papers ^3^2 H3 EI] an d their references). Their bulk density is given 

by 

£b(*)^(£^) +\ E (« 2 + J|0| 4 + fE(^) 2 + ^ 2 ' ( 4 ) 

a=l f3=m+l a=l 

where the position vector x has the representation (x 7 ) = (x Q ,x^) in Euclidean 
coordinates x 7 . We use the convention that a and /3 denote coordinate indices 7 
with 1 < a < m and m < (3 < d, respectively. Likewise, d a and dp mean the 
corresponding spatial derivatives <9 7 = d/dx^. At the level of Landau theory, the 
model has a continuous transition from a disordered to a homogeneous ordered phase 
for p > at r = 0. For negative p, a continuous transition to a modulated ordered 
phase occurs at a nonzero value of f. The transition lines between these phases merge 
at f = p = 0, which is an m-axial LP within Landau theory. 
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At LPs, systems exhibit scale invariance of the strong anisotropic kind: coordinate 
difference Ax a within the a-subspace scale as a nontrivial power (Axp) of the 
complementary ones Ax^, where 9, the anisotropy exponent, generally differs from 
1, and in Landau theory has the mean-field value MF = 1/2. 

Owing to this anisotropic scale invariance, boundary critical phenomena at LPs 
are richer than at critical points (CPs). For example, in the CP case, the orientation 
of the surface relative to the coordinate axes does not play any important role on the 
level of a description in terms of a <j> field theory. However, for systems at LPs the 
surface's orientation matters in an essential way, just as it does quite generally for 
systems with anisotropic scale invariance: Two fundamentally distinct orientations 
can be distinguished — one for which all a-directions are parallel to the surface (which 
we shall refer to as parallel), and another one for which the surface normal n is along 
one of the m a-directions (which we shall refer to as perpendicular). 

As a consequence of the different scaling behaviour of distances along a- and 
/3-directions, it depends on whether the orientation of the surface is parallel or 
perpendicular which boundary contributions C\ are potentially infrared relevant below 
the upper critical dimension d*(m) = 4 + m/2 and hence must be included in the 
action. The problem of constructing semi-infinite extensions of the models with bulk 
density J3J for d = d*(m) — e that are "minimal" in the sense that all irrelevant 
and marginal boundary contributions not compatible with the presumed 0(n) and 
Euclidean^ symmetries are discarded was considered for the case of parallel surface 
orientation in references |15lllfi| . where it was found that a contribution oc J] a (9 a ^) 2 
had to be included in the corresponding surface density L\ = c\, in addition to the 
one in equation J2J). The so-obtained semi-infinite model, defined by equation Q in 
conjunction with the boundary density 

o ? m 

4 = ^ 2 + ^E( 5 ^) 2 ' ( 5 ) 

was then utilized in references and |17j to determine the e expansion of the surface 
critical exponents of the corresponding ordinary and special transitions to second and 
first order, respectively. This extends or complements previous work based on the 
mean-field approximation |18| and Monte Carlo simulations for the axial-next-nearest- 
neighbour Ising (ANNNI) model |T5|. 

The model is a natural and simple-looking generalization of the semi-infinite n- 
vector model defined by equations (J2J and J2J, to which it reduces when m = 0. The 
fluctuating boundary conditions 0E] one obtains for it, 

d n 4> = (c - Xd a d a ) 4> , (6) 
where d n (= d z ) is the derivative along the inward normal n, are of the Robin type, 
generalizing the familiar ones of the CP (m = 0) case to nonzero values of A.§ As 
shown in references I17| , the new boundary term oc A has the effect of moving 
the fixed points onto which the ordinary, special, and extraordinary transitions are 
mapped to a nontrivial value A* = 0(e ) of the renormalized counterpart A of A. 

X We presume that bulk and boundary contributions breaking the rotational invariance in the a- 
subspace may be ignored. When m > 1, spatial anisotropics lead, in particular, to non-isotropic 
bulk contributions of the form T aia2 a 3 a l (d ai d a2 ( f>) ' d a3 d ai 4>, where T aia2 a 3 a 4 is a tensor whose 
symmetry is that of an m-cube or lower. The 0(e 2 ) results of reference 1141 for the bulk model 
indicate that such terms are relevant in the infrared for e > 0. 

§ Here and below we utilize the summation convention that pairs of the same a- and /3-indices are 
to be summed from a = 1 to m and j3 = m + 1 to d, respectively. 
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In the following we shall consider the case of perpendicular surface orientation. 
Previous work on it either was restricted to an investigation of the corresponding 
ordinary transition of the semi-infinite ANNNI model on the level of the mean-field 
approximation |2(J1 12 1| . or else employed Monte Carlo simulations to study both the 
ordinary and special transitions of this model [TT?] . 

Our aim is to construct and analyze an appropriate minimal semi-infinite model 
for general values of n, m and d. This turns out to be a greater challenge and more 
interesting than in the case of parallel surface orientation. Since the classical equation 
of motion of the order-parameter profile now involves a differential equation of fourth 
order in d z , two rather than one boundary condition are required (aside from two 
analogous ones at z = oo). Furthermore, being an x a coordinate, z = x m now scales 
naively as /i -1 / 2 where /i is an arbitrary momentum unit making fixp dimensionless. 
Compared to the case of parallel surface orientation (where z ~ fi , power counting 
gives more potentially relevant or marginal surface terms. Recently one of us has 
suggested an appropriate surface density and the associated boundary conditions 
22 . In section 2, we recall this density and derive the fluctuating boundary conditions 
it implies. Using the latter, we show that the existing additional surface contributions 
that are compatible with symmetry and short-rangedness of interactions, and not 
irrelevant according to power counting, are redundant. 

In section [3] we determine the free propagator subject to these boundary 
conditions, at the LP for general values of the surface interaction constants. Owing 
to its rather complicated form, explicit perturbative RG calculation for dimensions 
d = d*(m) — e are difficult to perform with it . In section 01 we show, following a 
suggestion made in reference [22], how this problem can be by-passed to some extent 
in the case of the ordinary transition. The asymptotic large-scale behaviour at this 
transition can be argued to be described by a theory with boundary conditions 



Hence one can work with the corresponding simplified free propagator at the price 
of having to deal with correlation functions involving besides the field <p a l so the 
boundary operator d^<p. From the anomalous dimension of the latter the required 
single independent surface critical exponent f3^ d of this transition can be inferred. 
Performing a one- loop RG analysis, we compute in section the e expansion of this 
exponent to 0(e) for general values of m. The result is used to estimate its value for 
the case m = n = 1 of the three-dimensional semi- infinite ANNNI model. Section 
contains concluding remarks. Finally, there is an appendix explaining how the required 
one-loop integral was calculated. 

While it is quite common to investigate field theories with Dirichlet boundary 
conditions (see [2] [23 and their references), we are not aware of previous studies 
of field theories satisfying the boundary conditions Q, barring familiar examples 
of hydrodynamic equations for fluids [21] • It therefore not unlikely that the work 
described below might also be of interest for other problems. 

2. Boundary density and fluctuating boundary conditions 

The boundary density suggested in reference [22] is 



<j> = d n (f) = . 



(7) 



m— 1 




j <p 2 + b cf>d n cf> + [y (d a 4>) 2 + f {d a <t>)d n d a <p] + -± {d n 4>f 



(8) 
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with / = 0. We have added the term ex / since it cannot be ruled out. Contributions 
breaking the symmetry among the m — 1 a-directions parallel to the surface have been 
excluded. That different values should be allowed for the coupling constants Aj^ and 
Ay should be obvious because the surface breaks the symmetry between a-directions 
parallel to it and the z-direction. Besides the monomials included in equation J5J|, 
there are a number of other surface operators one has to worry about, namely 

1 = <j>dfo, 2 = (d n 0)d> , Q s = 4>dl4>. (9) 

Before discussing these, let us first derive the fluctuating boundary conditions that 
apply to the model defined by equations (Q, and (JSJ). 

To this end, we compute the variation STL of the Hamiltonian. By integrations 
by parts, one gets 

d 



sn= / 5d> J^ + £flg £a 



8(8^) 



d n - 



(d n 5<t>) 



04, d{d n <t>) 1 ~' l d(d 2 n cj>) 



m-l Rr± 
a = l 



l d{d a <t>) 



m—1 



d{8n4>) d{dl4>) 



^ 9a d{dJ a <t>) 



(10) 



Equating the expression in curly brackets to zero gives us the classical field equation 



a(d a d a ) 2 - pd a d a - dpdp + f + 



= 



(11) 



Doing the same with the expressions in square brackets of the surface integrals 
J^6<p[. . .] and J^Sd n 4>[. . .] yields the boundary conditions 



\adl + {b- p)8 n + c± - [A,, + (/ - a)d r , 



771—1 

E 

a=l 



dl}<P 







(12) 



and 



-bdl + x ± d n + b-(f+a)22€}<f> = Q 



(13) 



respectively. 

These boundary conditions hold in Landau theory. Yet, they remain valid inside of 
averages, for the same reason that the classical equation (|llf) does so. To show this, one 
can make a shift <p - ► 4>+& in the functional integral defining the generating functional 
Z[J] oc J T>4> exp (—H + Jsjj J ■ <p) of multi-point correlation functions (<f). . . <f>). For 
a <& independent of <fi, the functional measure <fi does not change, and one arrives at 
the equation 



8H- 



J *) =0 



(14) 



where 4> must be substituted for 8<p in STi. Here (.) j indicates a normalized average 
in the presence of the source J, i.e., with the weight exp (— TL + f m J •</>). From 
this result (|14|l the validity of equations ^TJ-JTSJ inside of averages can be derived 
in a well-known fashion by setting the source to zero, either directly or after taking 
functional derivatives with respect to it, and exploiting the arbitrariness of <fr at and 
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away from the boundary. As usual, the source term yields extra contributions located 
at coinciding points to correlation functions generated by functional differentiation 
of equation 1)14(1, Corresponding fluctuation corrections would result from additional 
surface source terms such as J\ ■ 4>. 

It is now easy to understand why the surface operators given in equation do 
not have to be included in the Hamiltonian: These monomials involve d^4> an d d^<p. 
We can solve the boundary conditions (|13fl and (|12f> for these operators and substitute 
the solutions into the monomials. The expressions that result for the surface operators 
0i, 02 and 03 arc linear combinations of the surface operators retained in the surface 
density |JSJ.|| Thus the effect of the boundary operators 0i, 02 and 03 can be absorbed 
by a redefinition of the surface variables of the surface density (|SJ ■ That is, they are 
redundant and may be discarded. 



3. Free propagator at the Lifshitz point 

We now turn to the calculation of the free propagator G(x,x') at the LP. To this 
end we split the position vector into components parallel and perpendicular to the 
surface, writing x = (r,z). We choose z to be x m - The component r then 
involves the m — 1 a-components {x\ 1 . . . , x m -i) = r < and the al — m /3-coordinates 
{xp) — (aim+i, . . . ,Xd) = t>. For the wave- vector conjugate to r = (r<,r>) we use 
analogous conventions, writing p — (p<,p>). 

From equation l|ll|) one easily concludes that the partial Fourier transform 
G(p; z, z'), defined by 

G(x,x') = [ G(p; Z) z')jr< r -^ with / = / -^-^ , (15) 

Jp Jp jRd-l {2TT) a 

obeys the equation 

[a (p% - dlf +pl+p (pI -dl)+ t] G(p; z, z') = S(z - z') , (16) 

provided the order parameter profile (4>{x)) vanishes (as it does in the disordered 
phase). This equation must be solved subject to the boundary conditions 

[ad 3 z + [b+{f- a)p% - p] d z + c± + A|| ^}g(p; z = 0, z 1 ) = , (17) 
[-ad 2 z + \±d z +b+(f + &)p%\ G(p; z = 0, z') = , (18) 

and the requirement that the correct bulk propagator Gb(p,z — z') is obtained as 
z, z' — > oo at fixed z — z' . Furthermore, in order that the highest derivative, d^, 
produces the 5-function in equation l|16fl . both G(p;z,z') and its bulk counterpart 
must satisfy the jump condition 

[adlG{p-z,z%=l+l = l. (19) 

Further, G(jp\ z, z') must be symmetric under exchange of z and z' because G(x, x') is 
the inverse of a symmetric integral kernel A{x, x') associated with the Gaussian part 
d d x d d x' <p(x) A(x, x') <p(x')/2 of the Hamiltonian (including boundary terms). 
To simplify our analysis, we restrict ourselves to the LP, setting p = f — 0. It is 
also convenient to set the variable a (whose renormalized counterpart a changes under 

| In expressions such as • ™Z\ ®a 4> an< i d n <P ■ 2 a=l &a 4>> the derivatives d a evidently can be 
made to act as (— d a ) to the left, by means of integrations by parts. 
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RG transformations) temporarily to unity. The dependence on it can be re-introduced 
whenever needed by elementary dimensional considerations. For instance, for the free 
propagator these lead to the relation 

G( P< , P> ;z,z'|a) =«7- 1 / 4 G(a 1 / 4 p<,p > ;a- 1 / 4 z,a- 1 /V|l) . (20) 

The bulk propagator at the LP can be computed in a straightforward fashion. 
One obtains 

dk - ikz 
27 (fc 2 + p 2 -) 2 + pi 



G h (p,z) = I — ; (21) 



C -K+|Z| 



- , , -[k_cos(k_|z|) + k + sui(k_|z|)] , (22) 

4k+k_(k^ + K_) 

where we have written the roots of the denominator of the required Fourier integral 
as ±k_ ± i/c_|_, with 



K± = 7fVv /p4 < +p2 > ±p2 <- (23) 

In terms of the linearly independent solutions 

~ K+Z cos(k_z) for j = 1, 

, „~ K+Z sin(/c_z) for 7 = 2, , . 

Wj(p,z) = { K+z ( 1 \ ' (24) 

e + cos(k_2:) for j = 3, 

e K+z sin(K_z) for j = 4, 

of the homogeneous counterpart of equation Ijl6|) with p = f = and <r = 1, the free 
propagator of the semi-infinite system can be expressed as 

2 

G(p; z, z') = 6{z - z 1 ) W 3 {p, z) Vj(p, z') ,+(z~ z') (25) 



where 



4 

E 

k=l 



V J (p,z r )=yC jk W k ( P ,z l ). (26) 



The coefficients Cjk are chosen such that Vj—i t 2(p,z) satisfy the two boundary 
conditions (|17|l and l|18l) with p = and (7 = 1, that G(p; z, z') and its first and second 
derivatives with respect to z are continuous at z — z' , and the jump condition 11911 is 
fulfilled (again with a = 1). 

It is convenient to split off the known bulk propagator, writing 

G(p; z, z') = G b (p; z, z') + G s (p; z, z') . (27) 

Here G s , the part due to the surface, is (at least) four times differentiable in z and 
z 1 , and does not contribute to the jump of d?G at z = z', which originates entirely 
from Gb- Thus the jump condition l|19l) is taken care of. We can decompose Gb in the 
same way as we did for G in equations l|25l) and l|26|). with corresponding coefficients 
Cj^ . Evidently, only Gj^ +2 with j,f = 1,2 are nonzero and can be read off from 
equation (|21|) . but Gj^? = 0. Since G s must not have exponentially increasing parts, 
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four of the coefficients Cjk, namely Cjk — cjj^ with k — 3, 4, follow immediately from 
our result (|2"T|> for Gt>: 

C13 = C24 = -: — — ; — T\ = — —_ = — ^23 ■ (28) 

The remaining four coefficients must be determined from the two boundary 
conditions (|T7|) and ijT%|) for V\ and V2 . A straightforward calculation yields 

°n = ^ [B+ +pI + 84 + L- C] , 
Ci 2 = C21 = ^ [B_ + L - C] , 

C22 = [B K + pIk 2 _ -k 2 _L + C k ) , (29) 

where we introduced the short-hand notations 

K< = ° h + fP< > c p< = c_l + A||P< , 

B+ = b p< (b p< + 44) , B_ = 6 p< {b p< - 4k 2 .) , S K = 6 p< [6 P< (k 2 . + 2k 2 ) - j£] , 
C = c p< (A ± + 2k+) , C K = c p< [2kLk+ - A± (k 2 _ + 2k%)] , 

L — 2°Xj_k + (k 2 _ + k\) , D = -B++pl+L + C. (30) 

The resulting lengthy expression for the free propagator for general values of the 
surface interaction constants will not be used in the rest of the paper. The explicit 
form of the free propagator that we shall actually utilize is given in the next section, 
see equation lJ35]l. 

4. Asymptotic boundary conditions at the ordinary transition 

We now turn to the analysis of the ordinary transition. Let us begin by considering 
first the simpler case of the ordinary transition at a critical point (CP), and recall some 
essentials of its RG analysis. The critical behaviour at this transition is described by a 
fixed point at which the Dirichlet boundary condition <p\z=o = holds. One can choose 
this boundary condition from the outset (already for the bare theory) by setting the 
interaction constant c of the corresponding surface density J3J to the value c or d = 00 
01 12fil I27| . To obtain the behaviour of TV-point correlation functions involving 
fields <p close to the surface, one can employ the boundary operator expansion (BOE) 

<T>,z) a C a „4z) (d„0) re >) , (31) 

z— »0 

where </> ren and (<9 n 0) rcn denote renormalized operators. The required independent 
surface critical exponent /3° rd , which together with the usual bulk critical indices yields 
the other surface critical exponents of this CP ordinary transition via scaling laws, 
can be inferred from the scaling dimension A[d n 4>] = /3J rd /^ of d n <t> or the behaviour 

C dn ,(z) ~ z ^" d ^ (32) 

of the coefficient function Cq^. 

This CP ordinary transition must be recovered when p > for appropriate choices 
of the surface interaction constants. It is therefore reasonable to expect that the 
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Dirichlet boundary condition will prevail when p — > 0, i.e., at the fixed point describing 
the LP ordinary transition. To show this and to find the required second boundary 
condition, note that the surface variables with the largest /i-dimension are c± and b. 
At the trivial Gaussian fixed point, their dimensionless counterparts cj_ = /i~ 3 / 2 c_i_ 
and b = p~ 1 b transform under scale transformations p — > p,£ as 

c± ->c ± (£) =r 3/2 c_L and b -» &(£) = tr 1 6 . (33) 

Since the ordinary transition corresponds to the case where order is suppressed near 
the surface, we can take c± to be positive. Thus c± — > +oo in the large length-scale 
limit £ — > 0. Further, unless fine-tuned, b is not expected to vanish. Equation 1)33(1 tells 
us that b(t)/c(£) approaches zero as £ — > 0. The other coupling constants appearing 
in the boundary conditions 1(12|) . (|13l) . (|17|l and (|18f) with p = have /i-dimensions 
smaller than both cj_ and b. Hence the ratio of the associated running variables and 
either c± or b also becomes small as £ — > 0. We conclude that both 0| z= o and <9„0 
must become small in this limit in order for these boundary conditions to hold. In 
other words, the boundary condition J7J should apply in the large-length-scale limit. 
In the following we exploit this idea by choosing them from the outset for the bare 
theory. 

Let us denote by Goo(p; z, z') the free propagator that satisfies the boundary 
conditions 

G 00 (p;z,z') = d z G 00 {p;z,z') =0. (34) 

z=0 z=0 

It can be determined in a straightforward manner. One gets 

A / A , i\ A i a sin(K_z) sin(K_z') f , is 

Goo (p; «, 2 = Gb P; « - 2') - G b (p; 2 + z' — " " e ■ 35 

2k_ k+ 

Information about correlations near the surface can again be obtained via the 
BOE. Because of the boundary conditions Q, the leading operator contributing to 
the BOE of <fi now is <9 2 0; we have 

re >,z) « C a ^(z)(^r>) (36) 

instead of equation ((31(1 . The rcnormalized operator will be defined in the next section. 
Suffice it here to say that its scaling dimension 

A[9^]=/3(° 1 rd ' x) /^2 ) (37) 

where z^l2 is a bulk correlation-length exponent, gives us the surface critical exponent 
/CJlT' 3 '" 1 '' , and that by analogy with equation ((32(1 we have 

C 9 ^(z) ~ zfa-tir^VV"**) . (38) 
5. RG analysis of the ordinary transition 

Focusing our attention on the theory with the boundary conditions Q, we introduce 
the (N + M)-point cumulants 

/ iV M v cum 

G^; M L f) = ( J] ^. II %M (rj) ) ■ (39) 
\ 1=1 j=l I 

To regularize their ultraviolet (uv) singularities, we employ dimensional regularization. 
Apart from bulk uv singularities known from studies of bulk critical behaviour (121 1131 
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I14| , they have (primitive) uv surface singularities originating from the surface part of 
the free propagator, i.e. the last two terms on the right-hand side (rhs) of equation l|35|) . 
We use the reparametrization convention of references (141 1161 IT7| to absorb the uv 
bulk singularities, introducing renormalized quantities and renormalization factors via 

= Z* /2 ren , 6- = Z a a, ua- m ' A = F m , e p e Z u u, (40) 

T-T L p = p 2 Z T [r + A T p 2 ] , (p- plp)^ 1 / 2 = pZ p p, (41) 

with 

r(l + £ /2)r 2 (l-e/2)r(m/4) 
m ' e ~ (47r)( 8 +— 2e )/ 4 r(2 - e ) r(m/2) ' 1 j 

where we have explicitly indicated the dependence on a. Two-loop results for the 
renormalization factors Z^, Z a , Z p and Z T can be found in reference The 
renormalization function A T was computed to one- loop order in reference |16| : the 
result will not be needed in the following. 

In order to absorb the primitive uv surface singularities, additional counter-terms 
with support on the surface are needed. We introduce a renormalization factor Z 2 
and renormalized surface operator (d 2 4>) ren via 

d^cp = [Z a Z+] 1/2 {dl<t>T n = Zl' 2 dl4?™ . (43) 

Power counting shows that besides the surface counter-term resulting from Z 2 , G^ 0,2 -* 
requires an additive renormalization. Let us write the source term of the action needed 
to generate insertions of the bare surface operator d 2 cp as d r J 2 (r) ■ d 2 <p(r). 
Then the ^-dependent part of the renormalized action can be written as 

/ A d -\\Z- 1/2 J 2 - d 2 n cp^ + p 1 ' 2 a- 5 ' i S 2 J 2 ] . (44) 

J<8 1 J 

This requires some explanation. The counter-term involving Z 2 absorbs uv 
singularities ~ In A; as usual, the latter manifest themselves as poles at e = in 
the dimensional regularization scheme we prefer to employ. That such poles indeed 
occur, can be seen from the one-loop graph of G 1 - 1 ' 1 ' shown in figure^ Its closed loop 




Figure 1. One-loop graph of G( 1,:L '(p, 21). The stroke indicates the derivative d z . 
Crossed and open circles indicate points on and off the surface, respectively. Thus 
the crossed circle in conjunction with the double stroke on the right represents 
the surface operator d^(p- 

involves the distribution 

G o(x,x) =a- 1 l A p-V 2 z- i+2 'g e {m)F m ^ , (45) 

in which z = a^ 1 / 4 p 1 / 2 'z denotes the dimensionless counterpart of z, while g^im) is a 
number given by 

g e {m) = F~\ f G 00 (p; 1, 1|* = 1) . (46) 
J P 
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Since z > 0, the power z _4+2c corresponds to the generalized function i^ 4+2e of 
reference j^Hl- Upon substituting its well-known Laurent expansion [2H| (see also the 
appendix of reference [2]), we obtain 

G 0Q (x, x)/F m , e = a" 1 / 4 go(m) ^ + 0(e°) . (47) 

The quantity g e (m) is computed in the appendix for general values of m. Our 
subsequent results involve its value at e — 0, which is 

, , „ 3 f 24 (27r) 1 / 2 (2m-5)(m-l)r(m/2) 
g {m)=2-m-- 2-m + -r - — 



. + 2 8r[(m+l)/2] 

+ - 2 F 1 [2, (m - l)/2; (m + 2)/2; -1]) , (48) 
m J 

and reduces to 

flb(l) = -9 (49) 

in the uniaxial case m = 1. 

In accordance with its indicated /x-dimension, the counter-term cx S 2 diverges 
~ VA. Its is analogous to the surface counter-term ~ A in equation (3.135) 
of reference [2] needed in the CP case to renormalize the corresponding correlation 
functions (4> ai ■ ■ ■ 4>a,N®n4>h a • • • d n 4>b M ) in the cutoff-regularized theory. Just as 
Coo ~ A, the counter-term cx S 2 vanishes in the dimensionally regularized theory. 
We therefore do not consider it any further. 

Suppressing the tensorial indices ai,...,&M, we denote the renormalizcd 
counterparts of the cumulants l|39l) as 

G <f n M > = Z -(^+^ Z -M/2 G{N ,M) (5Q) 

A standard way of reasoning yields their RG equations 



V- a a N + M M 

^ + pp 9 p + — 2 — ^ + Y m 



G(f n M) =0, (51) 

p—u^a.p.r J 

where j3 p are bulk /3-functions in the notation of references [161 117) . Both (3 p as well 
as j3 T vanish at the LP and will not be needed in the rest of the paper. The functions 
(3 u (e,u) and r) a (u) = —f3 a /a are known to order it 3 from reference |13j . 

As a consequence of these RG equations, we have for the coefficient function 

CW(z)> 



C dU (z) = . (52) 



Solving this at the LP p = t = gives 

C 9 ^{Z) ~ (^- 1 /2)" 2 */(«) : (53) 

where tj 2 = rj 2 (u*) means the value of t] 2 at the infrared-stable fixed point u* = 0(e). 
The result can be combined with equation (|38|l to conclude that 

Pu d ' ±} = (d - m - 2 + ?7L2 + mO + rf 2 + 46)v L2 /2 

= {3 L + (46 + V * 2 )is L2 /2. (54) 
Here /3l = v\j2.\d — 2 + ?7L2 + — 1)] is the usual bulk magnetization exponent. 
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The scaling behaviour of other surface quantities can be derived along similar lines 
by exploiting the RG equations (|51|l in conjunction with the BOE i|36|) . In this manner 
one can determine how the associated surface critical exponents can be expressed in 
terms of four independent bulk critical exponents, such as 77L2, 8, v\ji and <p, and the 
surface magnetization index /3^ A . 

As an example, let us consider the surface susceptibility Xii(p) = 
G^ 2,0 \p;z = 0, z' — 0), the Fourier p-transform of the response of (</>| z= o) with respect 
to a surface magnetic field hi. To characterize the asymptotic low-momentum 
behaviour of this quantity at the LP, we introduce the surface exponent rj\\ by analogy 
with the CP case via 

(si„ g ),, (p^-W ior P> =0, 
X 1 (p) ~ < 1 (55) 

An ^ J p-^\^ p> m-i forp<=0, 

where "sing" means singular part. (At the ordinary transition considered here, xii 
does not diverge but approaches a finite value.) 

To identify m , we apply the BOE l|36|l to both external legs of & 2 '°\p, z, z'), and 
exploit the behaviour l)53|l of the coefficient function together with the scaling forms 
of G'- 2 ' ) and G*- ' 2 ) implied by their RG equations. This gives, for example, 

G^\ P ;z,z') z ^ oP> -^ +e [{zpi){z%)]^ + ' m , (56) 

which tells us that 

ryj ord ' ±) = 59 - 1 + r /L2 + rf 2 . (57) 

Other surface exponents, such as the surface susceptibility exponents 711 and 71 can be 
expressed in a similar manner in terms of P^i^'^ and bulk exponents. In particular, 
the usual scaling law 711 = ^2(1 — »7||) for the surface susceptibility exponent 711 is 
found to remain valid, as suspected by Binder and Frisch \21\ . 

Next, we turn to the computation of the still unknown RG function Z^. 
Using the result l|47|l . one can determine it to one-loop order in a straightforward 
fashion. Application of the distribution —S"'(z) to the external legs of & 0,2 \p) and 
& x,1 '{p; zi) gives 

- 6"'(z), [GqcX] 2 ) - 6 dlG 00 t I . (58) 
and 

(-5"'(z), G 00 (p; zi,z) [doo^l] (p; z, 0)) = 3 [g 00 K] (z x ,0) , (59) 

respectively. The implied poles of both Gren^(p) an d Grln\p; zi) get cancelled if we 
choose 

n 4- 9 11 

Z 2 = l- 90 (m) ~^- Te + 0(u 2 ) , (60) 

where we utilized the fact that Z§ = \ + 0(u 2 ). The derivative —ud u of this function's 
residuum at e = gives us the RG function 772 • Substituting for u its fixed-point value 
u* = 6e/(n + 8) + 0(e 2 ) then yields 

nl = 9o(m) u* + 0[K) 2 ] = go(m) + 0(e 2 ) . (61) 



( 
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The result can be combined with the known e expansions of the bulk exponents in 
equation (|54|) to obtain 



^T d,±) = K2/2)[4-e + % *]+0(e 2 ) 



1 + 



4(n + 8) 



■9m(0) 



0(e 2 ) 



(62) 



In the uniaxial scalar case m = n = 1 of the ANNNI model, this simplifies to 

-i -l 



/3r a) =i-^+o( £ 2 )= 



i 



24 1 



(63) 



Setting e = 3/2 in the first and second expressions on the rhs (i.e., in the "direct series" 
and the [0/1] Pade approximant) yields the d = 3 estimates P^^'^ — 0.31 and 0.59, 
respectively, which is to be compared with Pleimling's Monte Carlo result 0.62(1) |19|. 
Though these numbers are encouraging, our present knowledge of the series (|63ll to 
just first order in e clearly is insufficient to produce estimates that are competitive 
in accuracy with this Monte Carlo value. We therefore refrain from giving further 
extrapolated values for other surface exponents. Experience with the bulk case ^3] 
suggests that much better field-theoretic estimates should become possible once rfc is 
known to 0(e 2 ). In view of the simplifications entailed by the approach developed 
here, such a two-loop calculation should not be too difficult. 



6. Concluding remarks 

In this paper we have extended previous field-theoretic work on boundary critical 
behaviour at m-axial LPs EH El E3 by studying the critical behaviour at the 
ordinary transition of a semi- infinite system that is bounded by a surface perpendicular 
to an a-direction. This geometry was the one considered in the earliest investigations 
[2011^1 of boundary critical behaviour at LPs. However, to construct an appropriate 
minimal field-theory model and analyze it in a systematic fashion by means of 
modern RG method below the upper critical dimension turned out to be quite a 
challenge, mainly because significantly more potentially relevant (short-range) surface 
contributions must be considered than in the simpler case of where the surface is 
perpendicular to /3-direction. 

Taking up a suggestion made earlier |22j , we found a way to by-pass the enormous 
technical difficulties one is faced with when having to carry along many surface 
interaction constants. The basic idea is to choose the boundary conditions JTJl the 
theory is expected to satisfy on long length-scales at the ordinary transition from 
the outset, showing that they correspond to a fixed point of the RG. While we 
have not investigated deviations from this fixed point associated with modification 
of these boundary conditions (e.g., finite values of the surface couplings c± and b), 
our results are completely consistent with the physically reasonable expectation that 
these boundary conditions are associated with a stable fixed point of the RG. 

The benefit of our procedure is twofold: (i) Used in conjunction with the BOE, 
the RG equations we obtained enabled us to derive scaling laws relating the surface 
critical exponents at this transition to bulk exponents and a single independent surface 
index, such as (3^^''^. (ii) Owing to the gain in technical simplification, an explicit 

one-loop calculation could be performed to obtain the e expansion of and 
related surface exponents to first order. Of course, one cannot expect to get precise 
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numerical estimates for the critical exponents at d = 3 from such an 0(e) result. 
Nevertheless, the extrapolated values we obtained for /Sl™' 1 '"'"' in the case of the three- 
dimensional semi-infinite ANNNI model by direct extrapolation to d = 3 and via 
a [0/1] Pade approximant agree reasonably well with recent Monte Carlo results 
(albeit different numbers could be obtained if other scaling-law expressions for /?L° ld ' ± ' > 
were employed in conjunction with the best estimates of reference |13) for the required 
bulk exponents). More reliable field-theoretic estimates should be possible on the basis 
of 0(e) results. The required two-loop calculation appears to be quite manageable. 

A more difficult task is to perform an analogous two-loop RG analysis of the 
corresponding LP special transition. Since this requires the identification of the 
associated fixed point in the space of surface coupling constants, one cannot avoid 
retaining the dependence on these variables (see the analysis of the special transition 
for the case of parallel surface orientation |17) . for comparison). 

From a more general perspective, the present investigation is, to our knowledge, 
the first one dealing beyond the classical level with a field theory subject to the two 
boundary conditions (JJJ. Both our way of investigating this field theory as well as 
our result that this boundary condition correspond to a stable fixed point of the RG 
might have applications in other contexts. 
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Appendix A. Calculation of one-loop integral 

Here we compute the number g e (m) introduced in equation 1)46(1. 

Let us denote the bulk propagator in position space as G b (x — x'). Taking into 
account that Gb(0) vanishes at the LP in dimensional regularization, we obtain from 
equations l|45[) and (|46|l the result 

9e(m)F m!e = G 0( )(x,x\6- = l)\ z=1 = -G b (r = 0, 2z\l)\ z=1 + K d _ m K m _ x J e (m) .(A.l) 

The first term on the rhs follows from equations (7), (A5) and (A6) of reference |13j . 
It is given by 

G b (0, 211) = 2 - 1 -™ 7r -( 6 +™- 2£ )/ 4 „ F(2 ~ e) w N1 . (A.2) 
V ' 1 ' r[(m-2 + 2e)/4)] V ' 

In the second contribution we have split off the factors originating from the angular 
integrations, with 

K d = 2{Aix)- d l 2 /T(d/2) . (A.3) 
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In the remaining integral we make the changes of variables p > — > P = V>lv\ an d 
p< — > p = p</ V2 to obtain 



h{m) 



dp>pt 



d—m— 1 



with 



Jo Jo 



a + 1 a — m 

a = 2(i-4-TO = 4-2e, (3= — — , v = — - — 



We have introduced 



for which 



/ , \l/2 

( Vi + p 2 ± 1) , 



(A.4) 
(A.5) 

(A.6) 
(A.7) 



k^k + =P, k\-k 2 _=2, (k+ -i/c_) 2 = 2(l-iP). 

Upon integrating by parts with respect to p and using sin 2 x = (1 — cos2x)/2, wc 
arrive at 

J £ ( TO ) = — / dPP y 
" Jo 



/>OG 

/ dpp a sin(2pfc_) e - 2pfe + 
Jo 

I™ dpp a sm 2 (pk_)c- 2 P k + 
k- Jo 

/>oo /'OO 

/ dPP"" 1 *; 2 + / dPP^/W 
Jo Jo 

+ dPP"- 1 (Vl + P 2 + l)^ 2 



2^ 



with 



/•OO 

j(°> = / dpp a e - 2 ^+ = 2- 2 ^r(2/3) fc; 2 ^ , 

Jo 

J« = ^ dpp a sin(2pfc_) c- 2 P fe + = 2- 3f3 ^1 [(1 - iP)- fj - (1 + iP)" 
Jo 2i 



(A.8) 

(A.9) 
(A.10) 



and 



I {2) = I™ dpp a cos(2pk-)c- 2 P k + = 2~ 3fi ^^- [(1 - iP)-" + (1 + LP)-"] . (A.ll) 
Jo 2 

The remaining p- integrations can be conveniently performed by means of Maple. 

After some rearrangements the result becomes 



/ £ (m) = -2 2 - 5 r(^i)r(4-2 e ) 

+ (3-2e)2 1 - m / 2 - £ 



r(2-m/2-e) / m + 2e 

r(3/2- £ ) cos ^ 

r(l-m/4-e/2) 



' 1 7r ~ 



r(3/2 + m/4- e/2) 



^[ Maple 8, a product of Waterloo Maple Inc. 
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m — 2e 



20? ' 

to- 1 T(2 - m/2 - e)r(m/2) / m + 2e 

COS ( 7T ; 



to - 2 ?n + 2 N 
2 F 1 [2-e,——;——;-l 



2 r(2-e)r(m/2+ 1/2) 

( TO 1 TO 

x 2 Fi 2 e, :1 ;- 

2 1 V 2 '2' 2 ' 

For the uniaxial case to = 1, the required integration is simple, yielding 



(A.12) 



J e (l) = -y/2 f dpp 1 
Jo 



" e sin 2 



(A.13) 



= 2 e - 7/2 T(4 - 2e) [4 + 2 e cos(^e/2)] 
in conformity with equation (|A.12fl . 

Upon substituting the above expression l|A.12fl into equation (|A.lf) and using 
equation l|45|) together with the definition l|42|l of F m ^, one obtains the result for 
go{m) given in equation 
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